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The dynamics of ultracold neutral atoms subject to a non-Abelian gauge field is investigated. In 
particular we analyze in detail a simple experimental scheme to achieve a constant, but non-Abelian 
gauge field, and discuss in the frame of this gauge field the non-Abelian Aharanov-Bohm effect. In 
the last part of this paper, we discuss intrinsic non-Abelian effects in the dynamics of cold atomic 
wavepackets. 

PACS numbers: 03.75.Hh 
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I. INTRODUCTION 

Gauge potentials, and gauge theories in general, are 
crucial for the understanding of the fundamental forces 
between subatomic particles. The simplest example of 
gauge potentials is the vector potential in the theory of 
electromagnetism [l[ , which is an example of an Abelian 
gauge field. Non-Abelian situations, where the gauge 
potential is a matrix whose vector components do not 
commute, are surprisingly scarce in Nature. Candidates 
so far have mainly been restricted to molecular systems 
which are largely approachable only through spectro- 
scopic means. Other systems are liquid crystals which 
show the required non-Abelian symmetries [H, 0| • 

An elegant derivation and description of the emergence 
of non-Abelian gauge potentials has been presented by 
Wilczek and Zee |5(. It was shown by these authors 
that in the presence of a general adiabatic motion of a 
quantum system with degenerate states, gauge potentials 
will appear which are traditionally only encountered in 
high energy physics to describe the interactions between 
elementary particles. Ultracold atomic clouds are par- 
ticularly promising candidates for realising such scenar- 
ios, since the access to physical parameters is, from an 
experimental point of view, unprecedented. Extending 
the ideas of Wilczek and Zee, it was recently proposed 
that properly tailored laser beams coupled to degener- 
ated internal electronic states can be employed to induce 
non-Abelian gauge fields in cold-atom experiments [f|. 
Alternatively, a non-Abelian gauge potential can be con- 
structed in an optical lattice using laser assisted state 
sensitive tunneling ■ With the implementation of these 
proposals, ultracold atoms would offer a unique testbed 
for the analysis of non-trivial non-Abelian effects on the 
quantum dynamics of multicomponent wavepackets. 

In this paper we investigate the wave packet dynam- 
ics of a cloud of ultracold atoms in the presence of a 
non-Abelian gauge potential. In Sec. |TT] we discuss how 
this undoubtedly rather exotic scenario can be envisaged 
in a sample of cold atoms where the internal electronic 
energy levels are addressed by laser fields with a nontriv- 
ial spatial phase and intensity distribution. This setup 



opens up a number of new scenarios for ultracold gases, 
allowing for the study of non-Abelian atom optics, which 
naturally ties together optical and magnetic effects. Re- 
markably, as shown in Sec. [TH even very simple laser ar- 
rangements may induce non-trivial cold-atom dynamics. 
As a first example of this non-trivial dynamics, we discuss 
in Sec. IIIII a possible optical tweezer experiment includ- 
ing a non-Abelian flux, for which the population trans- 
fer between internal states crucially depends on the path 
taken (non-Abelian Aharonov-Bohm effect). This effect 
resembles indeed what one would expect from scattering 
protons onto a non-Abelian flux line where the proton 
can be transfered into a neutron [8| . The tweezer experi- 
ment discussed in Sec. IIIII iust involves the internal-state 
dynamics, without exploring the rich dynamics resulting 
from the interplay between external and internal degrees 
of freedom in non-Abelian gauge fields. Sec. IIVI is de- 
voted to the analysis of this interplay. In particular, we 
show that the dynamics of cold-atom wavepackets can be 
significantly affected by intrinsically non-Abelian effects, 
which are crucially dependent on the initial momentum 
distribution of the wavepacket. We consider in particu- 
lar the relevant examples of wavepacket propagation and 
wavepacket reflection at an atomic mirror. Finally, we 
conclude in Sec. W\ 
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FIG. 1: The tripod coupling scheme forms two degenerate 
dark states with a non-adiabatic coupling. The three laser 
beams fli i = 1, 2, 3 are arranged as two counter propagating 
beams (f2i and Q2) and one beam (^3) (of double intensity) 
in the perpendicular direction. 
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II. LASER-INDUCED NON ABELIAN GAUGE 
FIELDS 

In a recent paper Q it was shown that a non-Abelian 
gauge potential can be constructed in the presence of 
nontrivial light fields coupled to degenerate electronic 
states of cold atoms. For this we consider atoms with 
multiple internal states, see Fig. (1). For a fixed po- 
sition r the internal Hamiltonian Hq(t) including the 
laser interaction can be diagonalized to give a set of 
4 dressed states |x„(r)) with eigenvalues e„(r), where 
n = 1,2,3,4. The full quantum state of the atom de- 
scribing both internal and motional degrees of freedom 
can then be expanded in terms of the dressed states ac- 
cording to |$) = z~2 n =i *n( r ) lx™( r ))- If there are two 
degenerate dressed states and we can neglect the transi- 
tions to the other states we obtain a coupled two level 
system of the form 



1 

2m 



(-iW - A) 2 + V + $ 



(1) 



where the 2x2 potentials are given by 



Vn + (Xn(r)\V(r)\ Xm (r)), (2) 

A„, m = iH{x„(r)\Vx m (r)). (3) 

$n,m = |- [ <VXn|VXm}+5^Xn|VXfc}<Xfc|VXm> j 



k=l 



(4) 



The reduced 2x2 matrix A is sometimes referred to as 
the Berry connection and is related to a curvature (an 
effective "magnetic" field) B as 

Bi = -e m F kl , F kl = d k Ai - diA k - l -[A k , A t ] . (5) 

Note that the term h^iki[Ah, A{\ = (A x A), does not 
vanish in general, since the vector components of A do 
not necessarily commute. In fact this term reflects the 
non-Abelian character of the gauge potentials. 

To construct a scheme of laser-atom interactions that 
leads to a U(2) gauge potential we need degenerate (or 
nearly degenerate) dressed states. Such a condition 
is fulfilled e.g. for the tripod system shown in Fig. [1] 
A truly non-Abelian situation emerges if the matrices 
{A x , A y ,A z , $} do not commute. For this it is necessary 
that the off-diagonal element ift(xi(r)|V%2(r)) is non- 
zero. The Hamiltonian of the tripod system reads in 
interaction representation as [1, 0, [1(3] 



i(o 1 |o)(i| + n 2 |o)(2| + n 3 |o){3|) +h.c, (6) 



Parameterizing the Rabi-frequencies with angle and 
phase variables according tofli—fl sin 9 cos (p e lSl , fi 2 = 
fl sin 9 sin <p e lS2 = ft cos 9 e lS3 , where il = 



-\/|f2i| 2 + I O2 1 2 + |^ 3 | 2 , the adiabatic dark states read 

\D{) = sin0e iS31 |l) -cos(/>e iS32 |2), (7) 
\D 2 ) = cos6lcos0e iS31 |l) +cos6>sin(/)e iS32 |2) 

-sin0|3), (8) 

with Sij — Si — Sj . 

The gauge potential depends on the gradient of the 
dark states: 

An = ft (cos 2 </>V5 23 + sin 2 0VS"i 3 ) , 



A12 = ficosfl ^-sin(20)VSi 2 -iV0j , (9) 
A 22 = H cos 2 9 (cos 2 <t>VS 13 + sin 2 0V S^) , 



and 



$11 = |- sin 2 9 (- sin 2 (20)(VSi 2 ) 2 + (V0) 2 
2m \4 

$12 = sin 9 ( J sin(20) V5i2 - ivA (10) 
2m V 2 / 



sin(26»)(cos 2 </>VSi 3 + sin 2 0V5 23 ) - iV9 



$22 = ^-(~ sin 2 (29) (cos 2 cj)VS 13 + sin 2 0VS 23 ) 2 
2m \ 4 



{V9f 



This provides a remarkable versatility. Recent advances 
in shaping both the phase and the intensity of light beams 
makes it possible to choose practically any shape of the 
gauge potential provided the corresponding light field 
obeys Maxwell's equations. This is certainly the case in a 
two-dimensional geometry, but also in three dimensions 
light beams can be tailored [TJ U3] • In the Abelian case 
a nonzero effective magnetic field is obtained if there is a 
relative angular momentum between the two light beams 
and the intensity ratio is spatially dependent [HI, [H, EH • 
Surprisingly, the generation of a non-Abelian gauge 
field does not require any elaborate shaping of the three 
laser beams employed. This is indeed the case if we 
choose the configuration shown in Fig. [TJ Three plane- 
wave laser beams are used. Two lasers of equal intensity 
are counter-propagating in the x-direction with wave vec- 
tor k while the third one (of double intensity) propagates 
in the y-direction also with a wave vector n. With this 
arrangement, <f> = 7r/4, and 9 = tt/4 in the expressions 
above. The resulting vector potential is of the form: 



whereas 



V + $ = 



A = Hk 



Vi 



-e y e x /\/2 
e x /V% — ey/2 







(11) 



(12) 
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By choosing the laser detuning such that Vs — Vi = 
H 2 k 2 I Am we obtain a scalar potential proportional to the 
unit matrix, V + $ = Vil. Therefore the scalar potential 
can be safely neglected as far as the wavepacket dynamics 
is concerned. 



III. NON-ABELIAN AHARONOV-BOHM 
EFFECT 

In Ref. [7J] it was proposed that non-Abelian gauge 
fields created in lattices can be employed to construct 
non-Abelian atom interferometers. However, the read- 
out of any non-Abelian atom interferometer may be 
crucially handicapped by the non-trivial interplay be- 
tween external and internal degrees of freedom in the 
wavepacket dynamics of atoms in non-Abelian gauge 
fields (see Sec. IIV[) . However, this coupling between ex- 
ternal and internal dynamics may be prevented by con- 
sidering atoms trapped in mobile optical tweezers. If the 
tweezer potential is strong enough, the system may be 
investigated in the so-called single-mode approximation, 
in which both components share exactly the same center- 
of-mass wavepacket. As a consequence, the non-Abelian 
gauge field will just affect the internal dynamics of the 
atoms. In the following we envisage an experiment in 
which a cloud of ultra cold atoms is trapped by an opti- 
cal tweezer under the conditions discussed above. When 

I 



Fig. [3] shows the final population difference between the 
two dark states for both paths as a function of the path 
length ks. It becomes clear that the outcome of choos- 
ing the L or R path can be very different. We stress 
that this effect is not directly linked to the appearance 
of off-diagonal terms in the corresponding matrices of 
the vector potential, but rather it is inherently due to 
the non-Abelian character of the matrices A x and A y . 
This effect is remarkably similar to the scattering of pro- 
tons onto a non-Abelian flux line, where a conversion 
of the proton into a neutron is anticipated 0] . A more 
complete picture is obtained by defining the pseudo spin 



moving in the xy plane the atoms experience the gauge 
potential given by Eq. (TTTj) . We consider the case where 
the atoms are moved in the x and y direction (Fig. [2|) 
along two different paths: (clock-wise, L) from (0,0) to 
(0, s) and then from (0, s) to (s,s); (anti clock-wise, R) 
from (0,0) to (s,0) and then from (s,0) to (s,s). The 
initial state of the atom is assumed to be a linear super- 
position of both dark states: 

|¥(0)) - oos(rj)|Di) + fdnfaJI-Da) (13) 

where r\ is the mixing angle, and tp is a relative phase. 
The dynamics of the two level system obviously depends 
on the initial state, but more importantly, the final pop- 
ulations of the two dark states depend on which path is 
taken. After performing the clock-wise path the atoms 
are in the state 

= e tA * s / h e iA y s/h \y(0)) = cf|Di) + c£\D 2 ) (14) 

whereas after performing the anti clock-wise path we 
have: 

\*r) = e lA y s/h e lA * s/R \y(Q)) = cf \Di) + c$\D 2 ). (15) 

Using the vector potential given by Eq. (TTTj) a straight 
forward calculation yields 

(16) 
(17) 
(18) 
(19) 



a , L.R L,Rx 

S(Cj ,c 2 ' ) as 

S x = ^-.{cic 2 * - ci*c 2 ) (20) 

S y = i(c lC2 *+ Cl *c 2 ) (21) 

S z = \(\ci\ 2 -\c 2 \ 2 ). (22) 

With the pseudo spin representation we can follow the 
rotation of the spin vector as a function of position along 
the different paths. This is shown in Fig. [1] where the 
spin vector is seen to follow circular paths whose orien- 
tation changes when the direction of the atoms in real 
space changes. The role of the initial state is now im- 
mediately clear. Only a superposition between \D±) and 



cf = e- 4KS cos(-^=)cos(7 ? )+ie l ^-' ls / 2 )sin(-^=)sin(r;) 

cf = ie- lKS sin(-^)cos(?7) + e t ^- KS / 2 )cos(-^)sin(r/) 

cf = e _4KS (cos(-^=)cos(?7) + ie* v sin(-^L) sin^)) 
V 2 v 2 

c r. = e -»/2( isin (^) cos (^) + e ^ C os(^:) sin(r))). 
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FIG. 2: (Color online). The envisaged experiment. An optical 
tweezer moves the cloud of atoms along the left (L) path or the 
right (R) path. The final dark state population will depend 
on which path was taken. 
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FIG. 3: The difference in the populations in the two dark 
states depends on which path is taken. Fig. (a) shows the 
total difference jc^| 2 — | c^" | 2 as a function of the path length 
whereas in Fig. (b) we depict | cf* | 2 — c^] 2 . We assume 
as initial condition r\ = 7r/4, and ip = n/2. 

|Z?2) will result in a different final state of S as a func- 
tion of taking either the L or R path. Note in contrast to 
the previously considered laser-driven population trans- 
fer for tripod atoms [tj [l(| here the non-Abelian dynam- 
ics is due to the time-dependence of the phases of light 
fields "seen" by moving atoms rather than due to the 
time-dependence of the intensities of laser pulses. 

IV. WAVEPACKETS IN FREE SPACE 

The non-Abelian Aharonov-Bohm effect is a striking 
example where the internal dynamics of a two-level sys- 
tem is highly nontrivial. A question which is not of- 
ten addressed in the context of non-Abelian systems 




FIG. 4: The pseudo spin trajectories depend on the initial 
state and which path is taken: a) Left path with r\ — tv/8, ip = 
0, b) Right path with r\ = n/8,<p = 0, c) Left path with 
77 = 7r/4, ip — 7r/4, d) Right path with r\ — 7r/4, tp = 7r/4. The 
spheres in each figure indicate the initial state (I) and the final 
state (F). The black circle is always the path first embarked 
on. In all cases we have chosen ks = 34.5. This will cause 
the spin vector to traverse the circular paths several times in 
each plane. 

is the dynamics of a wave packet. This situation is 
clearly more complex compared to the previous non- 
Abelian Aharonov-Bohm scenario where an adiabatic 
motion with respect to center of mass excitations and 
shape oscillations was assumed. We now have to fully 
take into account the coupled internal and external de- 
grees of freedom. 

In the following we discuss the evolution of a cold 
atomic wavepacket in the presence of a non-Abelian 
gauge field A — {A x , A v , 0}. We consider that the atomic 
gas is sufficiently dilute, and hence in this paper we ne- 
glect the effects of the interatomic interactions. We re- 
strict ourselves to the case in which both matrices A x 
and A y are space-independent, fn order to simplify the 
discussion below, we consider Aj — hnMj, with j = x, y, 

where k has units of wavenumber, and Af? = 1. We 
assume as well that the scalar potential may be consid- 
ered as a multiple of the identity matrix (as discussed 
above). Removing unimportant global energy shifts the 
Hamiltonian for a free particle becomes 

H = -^V 2 i - (m£ + My-pl (23) 
2m m \ ox ay J 

The atomic wavepacket can be represented by a spinorial 
wavefunction of the form 

$(r,t) = ( dpj^'H&t). (24) 
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Thus, we have 

H*(r,t) = [ dpH p (p)$(p,ty^/\ (25) 



where 



Hp(p)=^-1 + — (M xPx + MyPy 

2m to V 



(26) 



Hence for any given p the equation of motion ih&(p, t) = 
H p {p)<5>{p,t) yields $[p,t) = eiq>[iHp(£)t/h}$(p,t = 0). 
This evolution can be analytically obtained after diago- 
nalizing the matrix H p (p) for every p. 

We are interested in comparing the wavepacket evolu- 
tion in the presence of Abelian and non-Abelian fields. 
If the fields are Abelian, i.e. [M x ,M y ] = 0, then we may 
find a common eigenbasis for both operators, in which 
Mj =diag{A^, Xj}. As a consequence, the eigenvectors 

£± of H p (p) are independent of p, and the total wave- 
function is at any time a linear combination of the form 
l>(f, t) = <&+(r, + $_(r, f)£_, where 



dpe- l ^e l P- f/h y±(p- r^,i : 



0), 
(27) 

withyf 1 = Kk(X£, Xf), and <f>± — *2 m i' +^ T 7^- Hence, 
the wavepacket evolution can be considered as an inde- 
pendent scalar evolution for the wavepackets in each com- 
ponent. In particular, it may be easily shown that the 
center of mass position of the wavepacket &±(r, t) grows 
linearly in time with a velocity ((p) +ff iz )/m. Hence the 
two wavepackets tend to separate during the time evolu- 
tion. 

The picture changes completely if [M x ,M y ] ^ 0. In 
this case the eigenvectors of H p (p) do depend on the mo- 
mentum p considered, and hence the time-evolution of 
the wavepacket depends in a non-trivial way on the mo- 
mentum distribution of the original wavepacket. We an- 
alyze in particular the center of mass (CM) motion of the 
wavepacket. The x-coordinate of the CM after a given 
time t is better calculated in the momentum represen- 
tation: (x) t = (ihd/d Px ) t = (e iflt ' h ihd/dp x e- lflt l h ) , 
where we have employed the Heisenberg picture. One 
can then easily obtain that: 



(x) t = (x) + —( Px )o + ( e'°ih— 



-iO 



(28) 



where O = (Kt/m)(M x p x + M y p y ). The last term in the 
previous equation leads to non-trivial effects, which are 
easily illustrated by considering the particular example 
NL = &~, M„ = o ~ : 



(x) t = (x) + — (p x ) + — {(c 2 a x + sca z ) 

TO TO 

9 „ Dill sin q , 

s a x — sea, so, 



sin 2q 
~2q~ 



sin 2q 
~2q~ 



29) 




FIG. 5: Evolution of the center of mass coordinate (x) in units 
of V2Ap/h, as a function of r = A P K,t/\/2m, for A x — hK,a x , 
and A y = h,Ka z , r\ — 7r/4, <p — 0. The dashed line is the 
function /(r) = r. For short times the nonlinear evolution of 
the center of mass becomes clear. 



where c = p x /p, s = p y /p, q = ntp/m, and p 2 = p 2 x +p 2 - 
Let us consider an initial Gaussian wavepacket 



*(r) = *(r 



cos^e 
sin rje 



iip/2 
-iip/2 



where ^(r) is a Gaussian centered in x = y = and with 
the Fourier- Transform ~ exp(— p 2 /Ap 2 ). Then: 



V2Ap 



-erfi(r) 



sin2?7COS(p, (30) 



where erfi is the imaginary error function, and r = ^^ t. 
Note, that contrary to the Abelian case, we have two in- 
herently non-Abelian effects. On one hand, the evolution 
of the center of mass motion is in general a non-trivial 
non-linear function of time. However, for r 3> 1 a lin- 

r is recovered, i.e. there is a 



ear behavior 



characteristical transient stage where an inherently non- 
Abelian-induced non-linear CM evolution occurs (see 
Fig. [5]). On the other hand, contrary to the Abelian 
(or scalar) evolution, the evolution of the CM motion 
depends on the initial width Ap of the momentum distri- 
bution. This effect can be traced back to the dependence 
of the eigenstates £± on p. 

A third effect can be observed if we consider a Gaussian 
wavepacket with an initial < p y >ot^ 0. In that case, if 
< p x >q— 0, one obtains: 



T 



sin 2q 
~2q~ 



- 1 



sin 2rj cos <p. (31) 



Hence the x-dynamics depends on the momentum distri- 
bution in the y-direction, contrary to the case of Abelian 
evolution. 
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FIG. 6: (Color online) Total density after t = 10{2m/hn 2 ), 
for (a) Ap = 0.2hn and (b) Ap = 0.6Hk. At t = 0, r\ = 0, 
93 = and (p) = 0. In the strong non-Abelian case the wave 
packet expands asymmetrically. In an Abelian situation with 
a radially symmetric effective magnetic field the expansion 
would be symmetric. 



Note that the details of the momentum distribution 
play a very important role in the wavepacket evolution in 
non-Abelian gauge fields. Obviously, if | (p) | 3> k the non- 



< 



K, 



Abelian effects become negligible. But even if \(p)\ 
an Abelian evolution is recovered if Ap -C K, i.e. the non- 
Abelian effects are clearer for wavepackets which at t = 
are localized in space with uncertainties < 1/k. The lat- 
ter effect may be explained, because if Ap <c hn then H p 
may be (to a good approximation) simultaneously diago- 
nalized for all relevant values of p in the distribution, and 
hence again two separated wavepackets as those for the 
Abelian evolution are recovered. In addition, it is impor- 
tant to realize that the particular evolution also depends 
on the initial spinor configuration of the wavepacket (al- 
though this dependence is not inherently non-Abelian 
since it also occurs in the Abelian evolution). 

Fig. [5] shows the results of our numerical simulations 
of the wavepacket evolution for the Gauge field discussed 



above. Note that contrary to the usual Abelian (or 
scalar) evolution, there is a stark difference in the evolu- 
tion of the shape of the wavepacket for different values 
of the momentum spreading Ap/ n. 

The non-Abelian character of the gauge field leads also 
to interesting effects in the reflection of atomic wavepack- 
ets. Ultra cold atomic wavepackets can be reflected 
at laser or magnetic mirrors [16|, M, O- For typical 
situations the reflection of the center of mass of the 
wavepacket can be considered as specular, i.e. the angle 
of reflection of the wavepacket with the normal vector of 
the mirror is exactly minus the angle of incidence of the 
original wavepacket. Mathematically, the reflection can 
be considered as the superposition (in absence of mir- 
ror) of the original wavepacket and an image wavepacket 
travelling with opposite momentum and with a dephase 
7r. For the case of wavepackets in non-Abelian gauge 
fields, the effect of the mirror cannot be mimicked by 
this image picture (since contrary to the scalar case, a 
sinusoidal solution is not an eigenstate of H p ). As a con- 
sequence the intuitive specular-reflection picture must be 
revised in the case of wavepackets in non-Abelian gauge 
fields, even for the cases discussed below, in which both 
internal components experience exactly the same mirror 
potential. 

Fig. [7] shows the reflection of the wavepacket for Ap = 
hn (i.e. for momentum spreadings for which, as discussed 
above, the non-Abelian effects are significant). It is clear 
from the figures that the non-Abelian dynamics after 
the reflection is certainly not trivial. Remarkably, the 
center-of-mass position does not show in general a spec- 
ular reflection. Fig. [5] shows the sum of the angle of inci- 
dence and that of reflection for different incident angles 
in the non-Abelian regime. For usual scalar (or Abelian) 
evolution, this sum equals zero. However, due to inher- 
ently non-Abelian effects, this sum is significantly differ- 
ent from zero. Moreover, contrary to the usual scalar (or 
Abelian) evolution, the angle of reflection crucially de- 
pends on the absolute value of the incoming momentum, 
and on the momentum spreading of the wavepacket. 



V. CONCLUSIONS 

Summarizing, gauge fields may be generated using ap- 
propriate laser arrangements with atoms with degenerate 
internal states. Using a very simple laser configuration, 
spatially homogeneous but non-Abelian vector potentials 
can be generated. In spite of this spatial homogeneity, 
the non-Abelian character of the vector potentials can 
lead to a surprisingly rich physics for the wavepacket dy- 
namics of ultra cold gases. On one hand, the free expan- 
sion dynamics of wavepackets crucially differs from what 
would be expected in scalar (or Abelian) cases. In the 
latter, the wavepacket center-of-mass follows a linear de- 
pendence in time. In the presence of non-Abelian fields, 
the wavepacket presents a non-linear time dependence 
during a transient time. In addition, and again contrary 
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bility of the corresponding Hamiltonian, the non-Abelian 
gauge fields introduce a dependence of the dynamics in 
different spatial directions. The wavepacket reflection off 
an atomic mirror is also significantly distorted by the 
non-Abelian gauge field. In particular, the reflection of 
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FIG. 7: (Color online). The reflection dynamics of a non- 
Abelian wave packet compared to a zero gauge field situation. 
The reflection takes place at x — — 7 where a steep potential 
is envisaged (gray area). The parameters were chosen to be 
K = 1, Ap = 1 and initial momentum po = — -^(x + y). 
a) The initial density distribution of the atomic cloud. The 
initial momentum kick is indicated by the arrow, b) The non- 
Abelian path of the center of mass, the inner (green) path, 
for the reflection is clearly different from the standard wave 
packet reflection with k = (red outer path), c) A snap shot 
of the wave packet at the time corresponding to the mirror 
image with respect to the i-axis. For k = the reflected angle 
is the same as the incident angle, d) A snap shot of the wave 
packet at the same time as in c). For the non-Abelian case 
the reflection dynamics is highly non-trivial where the center 
of mass path no longer is described by an incident angle equal 
to the reflected angle. 



to the scalar or Abelian case, the center-of-mass dynam- 
ics crucially depends on the momentum spreading of the 
wavepacket. Moreover, in spite of the apparent separa- 



FIG. 8: (Color online). The reflected angle plus the incident 
angle, 9 r + 9i, as a function of the incident angle 9i. The 
deviation from the standard case, 9 r + 8i = for a non- 
Abelian system is clearly seen. The parameters were chosen 
to be k = 1, Ap = 1 and initial momentum |po| = 8. 



the center-of-mass ceases in general to be specular, and 
the angle of reflection depends on the incoming veloc- 
ity and the initial momentum spreading, which is dif- 
ferent from the standard scalar case. The complex in- 
terplay between external and internal dynamics should 
make difficult the read-out of non-Abelian interferom- 
eters. However, an experiment performed with optical 
tweezers, may allow for the analysis of non-Abelian ef- 
fects in the internal dynamics of the atoms. In partic- 
ular, we have shown that such an arrangement can be 
employed for the analysis of the equivalent of the non- 
Abelian Aharanov-Bohm effect, where the final internal 
state of the atoms crucially depends on the particular 
path chosen. 
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